Fractal Weyl laws for quantum decay in generic dynamical systems 
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Weyl's law approximates the number of states in a quantum system by partitioning the ener- 
getically accessible phase-space volume into Planck cells. Here we show that typical resonances 
in generic open quantum systems follow a modified, fractal Weyl law, even though their classical 
dynamics is not globally chaotic but also contains domains of regular motion. Besides the obvious 
ramifications for quantum decay, this delivers detailed insight into quantum-to-classical correspon- 
dence, a phenomenon which is poorly understood for generic quantum-dynamical systems. 
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Phase-space rules provide powerful universal relations 
for classical and quantum systems. A time-honored ex- 
ample is Sabine's law, originally formulated in the con- 
text of room acoustics, which can be cast into the relation 
''"dwell = 4 V/ vA for the mean dwell time of a classical par- 
ticle escaping from a container, expressed in terms of the 
volume V of the container, the area A of the opening, 
and the particle's velocity v py]. In quantum mechanics, 
Weyl's law approximates the number J\f{E) of states with 
energy < E hy the number of Planck cells which 
fit into the accessible phase-space volume of the corre- 
sponding classical system (here h is Planck's constant 
and d is the number of dimensions) For a quantum 
particle in a container (e.g., an electron in a quantum 
dot), J\f{E) oc E^^'^ therefore follows a power law with a 
strictly determined exponent. 

In reality, quantum systems are open, and the eigen- 
states acquire a finite lifetime — the resulting resonance 
states constitute a fundamental concept across many 
fields of physics. However, there are reasons to believe 
that the synthesis of both mentioned phase-space rules in 
open quantum systems can modify Weyl's law. Evidence 
in this direction is provided by systems with globally 
chaotic classical dynamics, which exhibit a fractal Weyl 
law Af{E^r) oc CrE^^^'^ for the number of resonances 
with En < E and lifetime Tn > r (where the cut-off value 
T only enters the shape function c^-) [1, 0, & 0] • Remark- 
ably, the number dn in the exponent is not an integer; in- 
stead, it is given by the dimension of the strange repeller, 
which for a chaotic system is a fractal [7] . While this ob- 
servation already proofed suited to instigate a paradig- 
matic shift of the study of resonances in open quantum- 
chaotic systems, its direct practical consequences are nec- 
essarily limited: typical dynamical systems are not glob- 
ally chaotic, which has profound consequences on their 
quantum dynamics [8]. 

Here, we show that fractal Weyl laws indeed apply 
much more generally to generic dynamical systems, for 
which regular and chaotic motion coexists in a mixed 
phase space The key which reveals the fractal Weyl 
law is to restrict the resonance counting to a window of 
typical lifetimes {r < Tn < r' ^ where r <C Tdweii <C r'). 



We arrive at this conclusion by a combination of semi- 
classical arguments (based on a tailor-made phase-space 
representation of resonance wave functions which avoids 
problems with their mutual non-orthogonality) with nu- 
merical results for a paradigmatic quantun i- dy namical 
model system, the open kicked rotator 

We start our considerations with some general obser- 
vations about quantum-to-classical correspondence. A 
basic ingredient in the derivation of the ordinary Weyl 
law in closed systems is the mutual orthogonality of en- 
ergy eigenstates, which is guaranteed by the hermitic- 
ity of the Hamiltonian (equivalently, the unitarity of 
the time-evolution operator). Quantum-to-classical cor- 
respondence can then be exploited, e.g., by using a ba- 
sis of semiclassically localized states \x) which occupy 
Planck cells in phase space x = (q, p). Resonance wave 
functions, however, overlap with each other because open 
systems are necessarily represented by non-normal (nei- 
ther hermitian nor unitary) operators. In open quantum 
maps, for example, resonances are associated to the spec- 
trum of a truncated unitary matrix j\4 = QFQ, com- 
posed of the unitary time-evolution operator F of the 
closed system and the projector Q = onto the non- 
leaky part of Hilbert space (see Fig. [T] for an example 
of an open quantum map). Because of the truncation, 
the eigenvalues jUn = exp(— — 7^/2) of Ai lie inside 
the unit disk of the complex plane, ensuring that the de- 
cay rates 7^ = l/r^ are positive [lH. But since M is 
not normal, the associated eigenstates — which now de- 
scribe the resonance wave functions — are not orthogonal 
to each other. This circumstance, which is the root for 
the scarceness of analytical tools in non-normal problems, 
complicates the task of exploiting quantum-to-classical 
correspondence for the purpose of resonance counting. 

We circumvent this problem by applying standard 
phase space methods to an alternative spectral decom- 
position, the Schur decomposition [13] M. — UTU^ , 
which delivers the same eigenvalues (occupying the di- 
agonal of the upper triangular matrix T) but associates 
to them an orthogonal basis set U. For definiteness as- 
sume that all eigenvalues are ordered by their modulus, 
|Mi| < 1^2! < ImsI < ... < |mm|. The first r rows Ur of 
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U then form a complete basis for the r fastest decaying 
resonance states. Quantum-to-classical correspondence 
can now be exploited by considering the Husimi repre- 
sentation [13] 7ir{x) = X]m=i l(^l'^m)P of this subspacc. 
This provides insight into the regions in classical phase 
space which support the quickly decaying quantum reso- 
nances. An analogous construction can be based on the 
opposite ordering, |/ii| > |/i2| > l/^sl > • • • > |mm|, which 
focusses on the slowest decaying resonance wave func- 
tions. Because of non-orthogonality, the resulting 'fast' 
and 'slow' Husimi-Schur representations carry indepen- 
dent information on the resonance wave functions. 

Figure [T] demonstrates the viability of this method for 
a paradigm of quantum chaos, the kicked rotator [9>] with 
unitary time-evolution operator 

which classically reduces to the (symmetrized) standard 
map = p -\- /ccos(27rg + np) mod 1,^' = q -\- p/2 + 
p^ /2 mod 1. The matrix dimension M = determines 
the inverse effective Planck constant [l5| , while the kick- 
ing strength k determines the nonlinearity. For k = the 
classical system is integrable, while for /c > 7 the dynam- 
ics is globally chaotic (resonances in the chaotic variant 
have been studied in Refs. [3, Iiq|). 

In Figure m the kicking strength is set to /c = 2, for 
which the phase space of the closed system is mixed, 
as shown by trajectory segments (black dots) in panel 
(a). The color-coded areas superimposed on the phase 
space indicate the classical initial conditions for escape 
after 1, 2, 3, or 4 iterations when an opening is placed 
8it < q < 0.2 (rdweii = 5). Panel (b) shows the ordi- 
nary Husimi representation of quickly decaying resonance 
eigenfunctions with < \/in\ < 0.1 (M = 1280). Panels 
(c) and (d) show the 'fast' Husimi-Schur representation 
(0 < |/in| < 0.1, M = 160 and M = 1280), while panels 
(e) and (f) display the 'slow' Husimi-Schur representation 
(0.98 < \fin\ < 1; M = 160 and M = 1280) [16]. 

The ordinary Husimi representation demonstrates that 
the quickly decaying resonance eigenfunctions are all lo- 
calized in the region of escape after a single iteration. 
However, since these eigenfunctions strongly overlap, 
Planck-cell partitioning their phase-space support signif- 
icantly underestimates their number (predicting :^ 25 in- 
stead of r = 44 short-lived states for M = 160, and 200 
instead of r = 568 short-lived states for M = 1280). 
In contrast, the fast Husimi-Schur representation clearly 
maps out the classical escape zones and uncovers that the 
domain of quantum-to-classical correspondence increases 
with increasing M. The slow Husimi-Schur representa- 
tion maps out stable phase-space regions that are clas- 
sically decoupled from the opening by impenetrable dy- 
namical barriers; these regions do not significantly change 
as M increases. 

By construction, Planck-cell partitioning of the sup- 
port of the Husimi-Schur representations accurately es- 
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FIG. 1: (Color) Quantum versus classical escape in a kicked 
rotator (kicking strength k = 2.0) with opening at < 
q < 0.2. (a) Classical escape zones in phase space {q,p), 
color coded according to escape after 1 (red), 2 (blue), 3 
(green), or 4 (magenta) iterations. The black dots are tra- 
jectory segments in the closed system, whose phase space is 
a mixture of regular and chaotic motion, (b) Husimi repre- 
sentation of short-lived quantum resonances with decay time 
Tn < 1/(2 In 10) (decay factor \pn\ < 0.1), for Hilbert-space 
dimension M = 1280. The other panels show the 'fast' 
(panels c and d, \pn\ < 0.1) and 'slow' (panels e and f, 
I/in I > 0.98) Husimi-Schur representation, introduced in this 
work, for M = 160 (panels c and e) and M = 1280 (panels d 
and f). 

timates the underlying number of resonances. The fast 
Husimi representation therefore uncovers a proliferation 
of anomalously short-lived states driven by the emerging 
quantum-to-classical correspondence — among the total 
count of M resonances, the fraction of short-lived states 
increases as M increases. On the other hand, the slow 
Husimi representation shows that the fraction of anoma- 
lously long-living resonances supported by classically un- 
coupled regions remains fixed, which corresponds to the 
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FIG. 2: (Color online) Fractal Weyl law for resonances in the 
open kicked rotator with a mixed classical phase space (kick- 
ing strength k = 2.0, opening at < < 0.2). (a): Fraction 
P(/i) of resonances with decay factor |/in| > corresponding 
to a lifetime Tn > — 1/(2 In /x). (b) Fraction Ptyp of reso- 
nances in a band of typical lifetimes (0.1 < |/in| < 0.98), as 
a function of dimensionless system size M. In this double- 
logarithmic representation, the fractal Weyl law results in a 
linear dependence with a negative slope. 



ordinary Weyl law for this effectively closed-off part of 
the system. Crucially, we are led to conclude that the re- 
maining fraction of resonances with typical lifetime (cho- 
sen here to satisfy 0.1 < |/in| < 0.98) decreases as M 
increases, and therefore cannot follow an ordinary Weyl 
law. 

A detailed understanding of the resonance distribution 
can be obtained by counting the resonances in fixed life- 
time windows and comparing these counts for different 
values of M. We adopt probabilistic terminology and 
proceed in two steps. In the first step we determine the 
fraction P(/i) = prob( | /in | > of resonances with life- 
time exceeding a lower threshold r = — 1/(2 In /i). This 
defines a monotonously decreasing function interpolating 
between P(0) = 1 and P(l) = 0. In the second step, we 
extract the fraction of resonances within an interval of 
typical lifetimes (0.1 < < 0.98), which follows from 
Ptyp = prob(|/in| e [0.1,0.98]) = P(O.l) - P(0.98) (our 
results and conclusions do not depend on the chosen win- 
dow as long as it stays in the range of typical life times). 

Numerical results for the system with opening at < 
q < 0.2 are shown in Fig. O Panel (a) shows P(/i) for 
various values of M. As a function of /i, P decreases very 
sharply at the two extreme ends of the graph. At /i :^ 
we witness the influence of extremely short-lived reso- 
nances, while at /i 1 we observe the states which have 
very long lifetimes. Applicability of the ordinary Weyl 
law would entail that modulo small fluctuations, P(/i) is 
independent of M, since the uncertainty-limited resolu- 
tion of phase space increases uniformly when the Planck 
cell shrinks. The plot, however, shows that the body of 
the function P drops as M increases. This is due to the 
proliferation of the short-lived resonances (/i ^ 0), whose 
relative fraction among all resonances increases with in- 
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FIG. 3: (Color online) Same as Fig. [21 but with opening 
shifted to 0.2 < ^ < 0.4. 



creasing M, in agreement with the expanding domain of 
support of the fast Husimi-Schur representation. In the 
region of long-living states (/i 1), on the other hand, P 
does not depend significantly on M, in agreement with 
the observed M-independent support of long-living reso- 
nances in the slow Husimi-Schur representation. 

Complementing these trends for short and long life- 
times, the body of P becomes flatter as M increases. 
As shown in Fig. EJb), the fraction Ptyp of resonance 
with typical lifetime therefore decreases with increas- 
ing M. The linear fit in this double- logarithmic plot 
demonstrates that this trend closely follows a power law 
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fractal Weyl law, Myp MPtyp ex 

Figure [3] shows that the fractal Weyl law remains in- 
tact when the opening is shifted to 0.2 < q < 0.4, so that 
it couples to a larger part of the regular regions in phase 
space [see Fig. ID^a)]. In this case Ptyp oc M~^-^^, and 
therefore A/typ cx M^-^^. Compared to the situation in 
Fig. [21 the fraction of long-living states is now reduced, 
in keeping with the shrunken size of the classically un- 
coupled phase-space region. The fast Husimi-Schur rep- 
resentation [Fig. [41(b)] shows that the newly coupled parts 
support additional short-lived resonances. Consequently, 
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FIG. 4: (Color) (a) Classical escape zones and (b) fast Husimi- 
Schur representation of short-lived quantum resonances (M = 
1280) for the kicked rotator of Fig. [J with opening shifted to 
0.2 < g < 0.4. 
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the states with typical hfetime are stih associated to the 
chaotic regions. 

For globahy chaotic systems, the fractal Weyl law can 
be understood by associating the proliferation of short- 
lived resonance to quasi-deterministic decay following 
classical escape routes [4]. This quasi-deterministic de- 
cay requires classical-to-quantum correspondence, which 
is lost exponentially on a time scale given by the Ehren- 
fest time r^; = In M . The probability to reside within 
the system decays exponentially, too, and is governed by 
''dwell- The power law for the fractal Weyl law there- 
fore arises from the combination of two exponential laws, 
based on the relation exp(— r£;/rdweii) = M^^/^^^^^i^ for 
the part of phase space where classical-to-quantum cor- 
respondence does not apply; this region supports M ex 
^d-i/Ardweii resonances (delivering an accurate estimate 
for du if the opening is sufficiently small [7] ) . For generic 
dynamical systems, our results confirm the association 
of short-lived resonances to quasi-deterministic escape 
routes, while the resonances of typical lifetime are now 
associated to chaotic regions dominated by sticking mo- 
tion, where classical power-law decay oc rule in place 
of the exponential decay [13] • This can be reconciled with 
the fractal Weyl law when one assumes that the loss of 
quantum-to-classical correspondence in these regions is 
similarly modified into a power law oc t^, so that the 
Ehrenfest time takes the algebraic form te cx 
The fractal Weyl law then arises from the combination of 
two power laws, based on the relation r^" oc for 
the part of phase space where classical-to-quantum cor- 
respondence does not apply and dynamics is dominated 
by transient sticking motion. 

In summary, we have shown that typical resonances 
in open quantum systems obey a modified, fractal Weyl 
law, a phenomenon previously associated only to the 
non-generic case of systems with a globally chaotic clas- 
sical limit. We unravelled this law by introducing the 
concept of a Husimi-Schur representation, a phase-space 
representation of resonance wave functions which cap- 
tures maximal information on quantum-to-classical cor- 
respondence (circumventing the problem that resonance 
eigenfunctions are not orthogonal to each other). The 
formation of anomalously short-lived states that drives 
the departure from the ordinary Weyl law originates in 
quasi-deterministic decay along classical escape routes, 
whose phase-space support expands as one approaches 
the classical limit. The fractal Weyl law emerges from 
the interplay of the transient sticking of chaotic trajecto- 
ries to the stable components and the algebraic violation 
of quantum-to-classical correspondence, two mechanisms 
which are intimately related to a mixed phase space. 

The kicked rotator used here for illustration can be 
realized with atoms that are driven by pulsed optical 
waves [l9|. Experimentally, direct evidence of the frac- 
tal Weyl law is more likely to come from autonomous 



(non-driven) systems such as microwave resonators, in 
which advanced techniques allow for the accurate deter- 
mination of complex resonance frequencies [lO]; a mixed 
phase space is obtained for any generic smooth resonator 
shape. While the details of a mixed phase space consti- 
tute a unique fingerprint of a given dynamical systems, 
the general features (a hierarchy of stability islands em- 
bedded into chaotic domains where long-time transients 
arise from sticking motion) are remarkably robust. Based 
on this general phenomenology we expect that the frac- 
tal Weyl law for typical resonances is a generic feature of 
open quantum-dynamical systems. 
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